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We formalize the derivation of a generalized coarse-graining n-resolved master equation by intro- 
ducing a virtual detector counting the number of transferred charges in single-electron transport. 
Our approach enables the convenient inclusion of coherences and Lamb shift in counting statistics. 
As a Markovian example with Lindblad-type density matrices, we consider the Born-Markov-Secular 
(BMS) approximation which is a special case of the non-Markovian dynamical coarse graining (DCG) 
approach. For illustration we consider transport through two interacting levels that are either se- 
rially or parallelly coupled to two leads held at different chemical potentials. It is shown that the 
coherences can strongly influence the (frequency-dependent) transport cumulants: In the serial case 
the neglect of coherences would lead to unphysical currents through disconnected conductors. Inter- 
ference effects in the parallel setup can cause strong current suppression with giant Fano factors and 
telegraph-like distribution functions of transferred electrons, which is not found without coherences. 
We demonstrate that with finite coarse graining times coherences are automatically included and, 
consequently, the shortcomings of the BMS approximation are resolved. 

PACS numbers: 05.60.Gg, 73.23.-b, 72.10.Bg, 



I. INTRODUCTION 

The number of single electrons transferred stochas- 
tically through a small conductor in a given time in- 
terval obeys a statistics which is specific to the un- 
derlying transport process and to the details of the 
conductor— i^ii^iiii^i^. The deviations, especially for uncor- 
related charge transfer in single tunnel junctions (Pois- 
sonian statistics), from the Gaussian distribution, how- 
ever, are tiny and merely become visible in the tails of 
the distribution. Nevertheless, the exploration of the so- 
called full counting statistics (FCS) has established an 
active subfield of mesoscopic transport in recent years. 
Although theory work on FCS in mesoscopic transport is 
still highly dominating the field, several stimulating ex- 
periments were reported recentlyi'^'^ii^'^ii^ii^ii^ii^'i^iii. 
For example, the measurement of up to the 15th-order cu- 
mulant of tunneling through a single quantum dot (QD) 
and the observation of universal oscillations in FCS^^ in- 
dicate ongoing activities. 



Some recent theoretical work on counting statistics 

The theoretical study of FCS in quantum transport is 
mostly based on the computation of its cumulant gen- 
erating function. This turned out to be more conve- 
nient for practical purposes rather than the direct cal- 
culation of the distribution function. Various meth- 
ods were developed, e.g., the Levitov-Lesovik formula 
in the S-matrix formalismii^, diagonalization of Liou- 
villians of Markovian master equations^^ with the gen- 
eralization to frequency-dependent FCS by one of the 
authors in Ref.— , by means of Stochastic Path Integral 
formulation for classical stochastical network o^"'^^ , via 
a charge representation method^^, by an effective field- 



theorysii, with nonequilibrium- Greens functionssi, in a 
time-dependent Levitov-Lesovik approach^^, by a wave- 
packet formalisni^^, or for time-dependent FCS through 
the formulation by positive-operator-valued measure^l. 

The FCS for non-Markovian transport with gener- 
alized master equations were discussed in Ref.— and 
in Refi^ wherein the authors introduced an iterative 
scheme to compute the cumulants. The interrelation 
between waiting time distributions and FCS for single- 
particle transport were studied in Refi^S. Esposito et 
al^ have shown that the FCS for tunneling through tun- 
nel junctions obeys a fluctuation theorem relating distri- 
bution functions for forward and backward bias voltages. 

FCS and noise can be utilized as an important 
diagnostic probe of quantum coherence and decoherence 
mechanism a^^i^'^i'^^'^^i'^^i^''i'^^i'^^ . In principle, it turns out 
that the higher-order cumulants are very sensitive to 
coherent effects, as explicitly shown experimentally and 
theoretically for the second-order cumulant (noise) in 
Refia. 

The counting statistics of transport through single 
QDs with attached single phonon mode (nanoelectrome- 
chanical systems, Anderson-Holstein model) were ex- 
plored in RefSi ^-'^'^^i'^^'^^i^^'^^ and shows that the inelastic 
scattering processes strongly modify the current fluctua- 
tions. 

Further system-based theory on FCS can be found 
for bistable systems^ (generic FCS for telegraph cur- 
rent signals by means of stochastic-path-integrals), for 
Andreev scattering in an asymmetric chaotic cavity^, 
for a QD in the Kondo regime^, for Biittikers volt- 
age and dephasing probes by Pilgram et al^, for spin 
transfer through ultra-small quantum dots^, for trans- 
port through a molecular quantum dot magneli^, for 
chaotic cavities with many open channels^, for Andreev 
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reflection^, and for spin-transport with ferromagnetic 
leads^. 

Belzig studied the FCS of super-Poissonian electron 
transfer (bunching) in a single multilevel QD caused by 
dynamical channel blockade^. 

The combined counting statistics of transferred elec- 
trons and emitted photons in a single QD reveals 
crossovers between bunching and anti-bunching statistics 
for both species^i^. 



Master equations for single-electron transport 

Since we utilize a quantum master equation (ME) we 
will provide a brief survey of related previous work. 

A common starting point for a systematic perturba- 
tive description of single-electron transport through lo- 
calized states attached to leads (QDs, molecules, short 
carbon nanotubes, et cetera) is the non-Markovian ME 
with Born approximation (system-bath factorization) of 
the reduced density matrix. The assumption that the 
time evolution of the system state only depends on 
the present state leads to the Wangsness-Bloch Red- 
field approacbS2iffi'^>^. A further Markov approxima- 
tion yields a ME which is widely applied in the transport 
context e.g. i n63i64, 65^66, 67^68 _ well-known shortcoming 
of this approximation is that the positivity of the density 
matrix is not guaranteed, see e.g. Refj^ and references 
therein. To circumvent that, a secular approximation can 
be carried out and the ME acquires a Lindblad form^S. 
Another way to avoid the non-positivity is the so-called 
singular-coupling limit where the secular approximation 
is not necessary^2i2i. Here, we will present another alter- 
native based on (dynamical) coarse graining. 

A different access to the dynamics of the system den- 
sity matrix is provided by the Keldysh-contour formula- 
tion of Konig et al^^. This method allows for a system- 
atic diagrammatic expansion in the tunnel coupling. The 
noise for cotunneling though a single QD were studied in 
Refill. The FCS is obtained here by putting the counting 
fields (e*-^) in the tunnel Hamiltonian by hand^Si. 

Microscopic rate equations for transport through cou- 
pled QDs were derived in Refsi^i^^. These works are con- 
sidered to be benchmarks for the study of serially coupled 
QDs since the MEs turn out to be exact for infinite bias. 
For the first time, Gurvitz et alJ^ have been formulated 
the ME such that the system density matrix is resolved 
with respect to the number of transferred electrons. The 
FCS is then readily obtained. 

In RefiS Pederscn and Wackcr include broadening ef- 
fects in the ME description which become important for 
finite bias. They take into account the time evolution 
of coherences between different fc-states in the leads due 
to the tunneling processes. The numerical evaluation of 
the fc-resolved dynamics becomes rather expensive. How- 
ever, the results agree with state-of-the-art methods^^. 
For that technique FCS and noise is not considered yet. 

Recently, Leijnse and Wegewijsi^ reported on a ME 



approach for the reduced density using a Liouville-space 
perturbation theory. They systematically expand an ef- 
fective Liouvillian in Laplace space with respect to the 
tunnel coupling. The FCS in that framework has been 
studied in Ref,^. 



This work 



We address the single-electron transport through cou- 
pled QDs cither in parallel or in serial configuration 
for arbitrary Coulomb interaction strengths between the 
QDs and finite bias voltages. We use a ME based on 
lowest-order tunnel coupling. The approach is thus per- 
turbative and cannot capture Kondo effects. The dy- 
namical coarse graining (DCG) method^ prevents non- 
positive density matrices. For infinite coarse-graining 
times we recover the known Born-Markov-Secular (BMS) 
approximation whereas for small coarse-graining times 
the exact (non-Markovia n^^i^^ ) dynamics is obtained. 
By introducing a virtual detector for the transferred 
electrons at one tunnel junction we are able to calcu- 
late the time-resolved FCS and the noise spectrum at 
one junction. We will show that for both QD con- 
figurations the coherences in the energy-eigenbasis (off- 
diagonals of the system density matrix) play an impor- 
tant role and therefore cannot be neglected. Particu- 
larly, in the parallel setup interferences can lead to strong 
exponential current suppression with negative differen- 
tial conductanc e ''^'^'^1^^ . This accompanies with giant 
super-Poissonian Fano factors, Dicke-like noise spectra, 
and broad distribution functions of transferred electrons 
caused by the telegraph-like current signal. For the 
serial setup the BMS approximation (vanishing coher- 
ences in the energy eigenbasis) leads to unphysical behav- 
ior when the tunnel coupling between the dots is small 
and can be fixed by assuming a finite coarse-graining 
time. In our computations the Lamb shift terms (level 
renormalizationSS) arc always included and their effect is 
visible when populations and coherences in the energy- 
eigenbasis couple. 

In Sec. |TT] we show how we introduce a virtual detec- 
tor to the QD system and how we obtain the FCS and 
the noise spectrum. The coarse-graining method and the 
BMS approximation are presented in Sec. IIIII The mod- 
els, two QDs in series and in parallel are introduced in 
Sees. IIV Al and IIVBI respectively. Sec. |V] contains the 
equilibrium results, Sec.|Vl]the results for the serial con- 
figuration and Sec. IVIII the results for the parallel setup. 
In the Appendices \K\ and IB] we provide the DCG Liouvil- 
lians for parallel and serial setups, respectively. The eval- 
uation of the Lamb shift terms appears in Appendix [Cl 



3 



II. COUNTING STATISTICS BY A VIRTUAL 
DETECTOR 

When one wants to describe counting statistics with 
a master equation approach, one is first faced with the 
problem that the number of tunneled particles is actu- 
ally a bath and not a system observable. In addition, for 
systems with coherences it appears not so trivial to iden- 
tify matrix elements of the Liouville superoperator with 
jump superoperators. Therefore, we perform the count- 
ing statistics by adding a virtual detector with infinitely 
many eigenstates in e.g., the right lead. Formally, this is 
done by modifying the tunnel Hamiltonian^ 

di (g) 4^ ^ (g) 6^ (g) cj.^ , (g) cjj^ ^ di (g 1 (g) 4^ , (1) 

where di/dl and Cfea/c^^ denote the annihila- 
tion/creation operators in the dot system and the 
lead a = L, R, respectively. The detector excitation 
operator 

-t-oo 

6t= ^ \n + l}{n\ (2) 

n—~oo 

increases the occupation of the detector by one each time 
an electron is created in the right lead. When we for- 
mally treat the detector operators b and as system op- 
erators, we may perform the conventional Born-Markov- 
Secular approximation as presented e.g. in Refi^ and 
obtain a Lindblad^^ type master equation. For example, 
assuming a two-level dot system characterized by Fock 
states \n1n2) {rii G {0, 1}) the Lindblad form would op- 
erate in an infinite-dimensional Hilbert space spanned by 
|00)g)|n), |01)(g|n), |10)(g|n), and |ll)(g|n). We can gen- 
erally decompose the density matrix in this Hilbert-space 
as 

+00 

ps{t)^ ^ p("™)(i)®|n)(m| , (3) 

n,m— — 00 

where act on the dot Hilbert space spanned 

by |00), |01), |10), and can thus be represented 

as 4 X 4 matrices. By taking the "matrix elements" 
p^"\t) = — {n\ ps{t) \n) we see that the resulting 

n-resolved density matrices are related vic^ 

= LoP^"^ + L+p^"-^') + , (4) 

which implies that we can ignore dynamics of 
for m 7^ n. Note that the above approach also captures 
the small bias range leading to bidirectional transport, 
which leads to the occurrence of both and p'"^^' 

in contrast to other commonly used n-resolved master 
equations22i^i2Si2S. The generalized Liouvillian superop- 
erator Lq contains transitions between the system and 
the left lead, whereas L+ corresponds to jumps from the 
system towards the right lead, and vice versa for L_. 



Now, according to the measurement postulateSi count- 
ing m particles at time t in the detector would project 
the density matrix ^ to 

such that we may interpret 

P„W^Tr{/")(t)} (6) 

as the probability that n particles will be found in the 
detector when we measure at time t. Note however that 
the n-resolved density matrices p*-"-* (t) derived in this way 
are still positive semidefinite, since Eqn. ^ corresponds 
to a Lindblad form in a higher-dimensional Hilbert space. 

The n-resolved master Eqn. ([4]) can be Fourier- 
transformed pix^t) = P''"^^)^™'^' such that we can 
consider the Markovian cumulant generating function 
SMKix^t) for the probability distribution ^ 

gSMK(x,t) ^ J2 Pn{t)e"''' = Tr je^^'^^V} , (7) 

n 

where L{x) = Lq + e+''^L+ -I- e^^^L- and the initial 
density matrix is conventionally chosen as the station- 
ary state (fulfilling L{0)p = 0). If there is a distinct 
eigenvalue of the Fourier-transformed Liouvillian with a 
largest real part Ao(x)j this eigenvalue will dominate the 
evolution for large times t 

gSMK(x.t) ^ Tr{pOe^°(x)*+pie^i(x)* + ...} 

« e^o(x)tTr{pO} =e^«(^)*, (8) 

where Ao(0) = (stationary state). In this limit, the 
derivatives of the cumulant generating function can be 
approximately determined from 

i-id^)'SMKix,t)\^^,^ {-id^)'Xoix)\^^ot. (9) 

For the Markovian stationary current we obtain from 
the cumulant generating function 

/=-zTr{L'(0)p} , (10) 

where conventionally p denotes the (normalized) station- 
ary state of L{0), see also22i2^ for similar expressions. 
The Markovian finite-frequency noise can be obtained 
from the second cumulant via the MacDonald noise 
formula^! 

00 

Sniuj) = J wsin(c.t)^ (^{n'it)) - {n{t)f) dt , (11) 


where the regularization 

uj sm{ujt) hm [a; sin(a;t) -|- e cos(wi)] e"*^* (12) 



4 



is implied^^. It is straightforward to show that the sec- 
ond term in the McDonald formula pTjl is given by 
{n{t)) = It with the Markovian current pO|) and one 
can analytically obtain the associated integral. The first 
term in Eqn. (|lip can be written as a Laplace transform, 
such that one finally obtains for the frequency-dependent 
noise 



2L'(0) 
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iujl + L{0) 



L'(0)-L"(0) 



(13) 



where L'(0) ^ d^L{x\^, and L"(0) ^ d^Hx)]^^^, 

see alscM^ for similar relations. Note that L{0) is sin- 
gular, such that in order to obtain the zero-frequency 
limit of the above expression, one may either evaluate 
the pseudo-inverse of L{0) or directly deduce it from the 
dominant eigenvalue 5(0) = — Aq(0). 



III. COARSE-GRAINING AND THE BMS 
APPROXIMATION 

A. Coarse-Graining Master Equation 

It is well known that for nontrivial systems the conven- 
tional Born-Markov approximation scheme does not nec- 
essarily lead to Lindblad-type master equations. In or- 
der to obtain these, an additional secular approximation 
(termed BMS throughout this paper) is required. Alter- 
natively, the singular coupling limit also yields Lindblad- 
type master equations^^iZL. 

However, it is also known that Lindblad type master 
equations can also be obtained using coarse-grained time- 
derivatives2Si21. The Liouvillian will then depend on the 
coarse-graining timescale. For example, one may match 
the perturbative solution to the equation 



dt 



p{T,t) = L{T)p{T,t) 



(14) 



for the reduced density matrix (using that L{t) is small 
in the interaction picture) with the perturbative second- 
order solution of the von-Neumann equation for the 
full density matrix in the interaction picture at coarse- 
graining time T 



s^M>g = TrB{C/(r)pg®p^C/t(r)} 



(15) 



In the above equation, U{t) denotes the time evolution 
operator in the interaction picture (including time order- 
ing). For a decomposition of the interaction Hamiltonian 



M 



(16) 



fines the corresponding Liouville superoperator as 



L{r)[p] = E - / '^^i^* 



-C^^{tut2)eit2-ti)pAlih)Al{t2) 



+ Ca0{tut2)Ai3{t2)pAl{ti 



(17) 



which can be shown to be in Lindblad forni^^. In the 
above equation, the time-dependence arises from the in- 
teraction picture and the generalized bath correlation 
functions have been introduced as 

Ca(3{tut2) = TTB{Ba,{ti)B,3{t2)pB} , 

Csp{ti,t2) = Tib {Bl{ti)Bf3{t2)pB} , 
C^0{h,t2) = Tre {Bl{h)Bl{t2)pB} ■ (18) 

With the relation S{x) = 5 [1 -|- sign(a;)] one may sep- 
arate the Liouville superoperator (jl7p into dissipative 
terms of Lindblad form^S and non-dissipative terms 
(Lamb-shift) that can be expressed by a commutator 
Li^s{t)p ~ — i [iJcff (t), p] with an effective Lamb-shift 
Hamiltonian Hcs{t) = iJ^jj(T). For equilibrium baths 
{[Hb,Pb] = 0) the bath correlation functions (fT8|) will 
only depend on the difference of their time arguments and 
we may introduce the even (7ij) and odd (cry ) Fourier 
transforms 
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+00 



Cy(t)sgn(i) = ^ y" a,jiLo)e+'^'dLO (19) 

—00 

of the bath correlation functions. With these, the two 
time integrations in Eqn. (jl7p may be performed ana- 
lytically and just the one-dimensional integral over lu re- 
mains. 



B. Coarse-Graining Schemes 

In order to go beyond conventional Markovian ap- 
proaches, the coarse-graining time r must scale with the 
physical timo^. For example, in order to approximate 
the short-time dynamics of the exact solution well (and 
thereby also the dynamics of the non-Markovian mas- 
ter equation), r should scale linearly with the physi- 
cal time. In contrast to conventional Markovian mas- 
ter equations, the dynamical coarse-graining method 
(DCG) ^^i^^ therefore formally solves all coarse-graining 
master equations ()14p and then fixes the coarse-graining 
time parameter r = t in the solution 



into system (Aa) and bath (Ba) operators Eqn. de- 



Pit) 



(20) 
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For the pure-dcphasing spin-boson model, this choice 
yields the exact solution for the reduced density matrix. 

In addition, one can show analytically, that in the 
limit T —^ oo (an equilibrium bath assumed), the Born- 
Markov-secular approximation is reproduced (see^ for 
the detailed proof): The result of the temporal inte- 
grations in (jl7|) may be phrased in terms of band filter 
functions sinc(a;) = sin(a;)/x, which converge to Dirac 5 
distributions for r ^ oo, see appendix [XI for an exam- 
ple. Thus, in this limit all integrals collapse and pT]) 
finally reduces to the BMS Liouvillian. In the short-time 
limit, we automatically approach the exact solution (and 
of course also the non-Markovian master equation) by 
construction and for all times positivity is automatically 
preserved. 

However, in the long-time limit it is not a priori clear 
why the coarse-graining time r should always scale lin- 
early with the physical time t. In that case, we would 
always obtain the BMS stationary state in the long-time 
limit. This state has nice classical properties and yields 
the exact (non-perturbative) solution for the spin-boson 
pure dephasing model, but conflicts with some exact 
quantum solutions - as already exemplified by the single- 
resonant level mode l^^'^^ . Later-on, we will demonstrate 
further shortcomings of the BMS approximation. Es- 
pecially when one goes beyond the weak-coupling limit, 
a finite coarse-graining (FCG) time may lead to better 
results, such that one might alternatively investigate de- 
pendencies of the form 



r(t) = 



t 



1 



(21) 



Therefore, when one is only interested in stationary 
(Markovian) results, one may simply evaluate L(rniax)- 



C. Application to Transport 

It is rather straightforward to apply the coarse- 
graining method to rt-resolved master equations: 
Eqn. ([U will appear as 

p(")(T,t) ^ Lo(r)p("^(r,t)+L+(r)p("-i)(r,t) 

+L_(r)p("+i)(r,t), (22) 

and we can apply the same arguments that lead to 
Eqn. ([7]) to define a non-Markovian cumulant generat- 
ing function 



(23) 



where consistency with the BMS approximation is 
achieved by L{0) = lim L(0, r) and we have L{Q)p ~ 0. 

r — >oo 

When observables (cumulants) should be calculated, the 
actual scaling of the coarse-graining time with the phys- 
ical time is of importance, since from the above equation 
we can deduce the time-dependent cumulants by taking 



derivatives with respect to x ^-nd t. For example, we will 
refer to 



I{t) EE ^(r.(r(t),t)) 

= ^ (-z9^)W(x,T(t),t)|^^o (24) 

as the time-dependent current. When r is fixed (Marko- 
vian case), this current becomes time- independent and 
when T oo, we obtain the BMS current. 



D. Decomposing the Fock space 

It should be noted that at first sight fermionic hop- 
ping terms appear incompatible with the tensor-product 
representation of the interaction Hamiltonian (|16p : Op- 
erators acting on different Hilbert spaces commute by 
construction, whereas in the fermionic hopping terms the 
dot and lead operators anticommute (this feature is also 
exploited in the exact solution a^^i^^ ). However, we can 
map such hopping terms to a tensor-product representa- 
tion. From the original anticommuting fermionic opera- 
tors (denoted with overbars), such a representation can 
be obtained from the decomposition di = ^ 1 ® 1, 
3.2 = crl ® ^2 ® 1, and Cka = erf (g) erf (8) Cka, where 
= ^ (c^ -t- i(J^), which is similiar to a Jordan- Wigner 
transform. In this decomposition, the first two Hilbert 
spaces refer to the two-dimensional spaces of the first 
and second system site, respectively, and the last Hilbert 
space is simply the (infinite-dimensional) Fock space of 
the leads, within which the Cka operators obey the usual 
fermionic anticommutation relations. For the double dot 
systems considered here, the Hilbert space of the system 
is defined as the Fock space with up to two particles - 
spanned by the basis |00) , |01) , |10) , |11). In this system 
Hilbert space we can now define new fermionic operators 
via c?i = — (Tj^(3)(t| and ^2 = — l®cr^. The generalization 
of these mappings to systems with more than two sites 
is straightforward. Such mappings appear to be usually 
performed tacitly in the literature, see e.gj ^^'^^i^^ . 



IV. MODEL HAMILTONIANS 
A. Two Levels in series 
We consider the Hamiltonian^ii^i^iSl 



€id\di + e2<^2'^2 + Ud\did\d2 



H — Hs + HsB + H] 
Hs 

+Tc{didl + d2d\ 

HsB 



k 

\^kLdi ® 



,(25) 
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where d|^2 create an electron with different quantum action Hamihonian 



numbers on the dot, respectively, and cj,^y^ create elec- 
trons with momentum k in the left /right lead. The pa- 
rameters 6; denote the single-particle energies, U mod- 
els the Coulomb interaction, and Tc denotes the interdot 
tunneling rate. This spinless model could be motivated 
by a large magnetic field that leads to complete spin po- 
larization in the leads, such that only one spin would 
need to be considered, or - alternatively - by orbitals 
where all tunneling processes are completely symmetric 
in the electronic spin, such that it may be omitted from 
our considerations. We assume the symmetric bias case, 
such that the chemical potential for the left and right 
leads are kept at /^^ = +V/2 and fi^ = —V/2, respec- 
tively, where V denotes the bias voltage. The model 
is depicted in Fig. [T] The Hilbert space of the system 
can either be spanned by the localized (particle) basis 
|00) , |01) , |10) , or the hybridized eigenbasis of the 
system Hamiltonian 



A2 
A3 
A4 



di (g) 1 , 



k 
fc 
k 

= ^ tljiCkR ■ 



(27) 



The resulting expressions for the bath correlation func- 
tions and the Liouville superoperator become rather 
lengthy in the general case, but their derivation is out- 
lined in appendix |^ We assume a continuum of bath 
modes, such that we may introduce the tunneling rates 

Tl{uj) = 2^^|tfci|'j(c^-efcL), 

k 

TRiio) = 27:J2\tkR\''S{Lu-ekB.), (28) 





- 100) , 




1 - 


E-) 






^/n - 




1 - 


E+) 




- 


\E3) 


- Ill) , 



5 + v/JM"47f) |01) + 2T, 1 10) 



-2r,|01) + (S+^S^+AT^) |10) 



(26) 



with the normahzation N ^ 4T'^ + [5 + ^PTlT^j 

and the single-particle splitting (5 = ei — £2- It can be 
shown easily that for energetic degeneracy (5 = the 
energy eigenbasis remains non-local even in the limit 
Tc — *■ 0. We will see that in the BMS approximation this 
leads to unphysical artifacts for some parameter values. 

With introducing a virtual detector as described in 
Sec. ini we identify four coupling operators in the inter- 



which we will assume to be approximately constant (flat 
band limit). Note however, that for the evaluation of the 
Lamb-shift terms it is necessary to assume a cutoff, see 
the discussion in appendix O 

Generally, we observe that for finite coarse-graining 
times T, the Liouville superoperator (jl7p couples the 
six matrix elements poo,oo, P- , Pii.ii, 
to each other. In contrast, in the BMS limit (t 00) 
and assuming a nondegenerate spectrum, the popula- 
tions poo.oo, P-.-, /0+,+, Pii,ii hi the energy eigenbasis 
will decouple from the coherences such that 

it suffices to consider a 4 x 4 Liouville superoperator. 
This is a general property of the secular approximation^ 
(also rotating wave approximation^!) . When the single- 
charged states become energetically degenerate ei = £2 
and Tc = 0, this decoupling does not take place in the 
BMS limit. 



Two Levels in Parallel 




FIG. 1: [Color Online] Depending on the bias voltage V, par- 
ticles may travel from left to right or vice versa through the 
effective double dot system. Setting one tunneling rate {ttL, 
tkR, or Tc) to zero inhibits the current. 



e.g 



We consider the Hamiltonian (compare also 

36. 37.71. 83. 93^ 



H = 
Hs = 
Hb = 

HsB = 



eid\di + e2d\d2 + Ud\did2d2 , 



[^kLd\ ® CkL + tl^dl ® CkL 

k 

® CkR + tlj^dl ® CkR + h.c. , (29) 



where we now have four tunneling rates instead of two for 
the serial model As before, we omit the spin from 
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our considerations. Note that in this model, there is no 
interdot hopping, but the model becomes nontrivial due 
to the Coulomb interaction term. The model is depicted 
in Fig. H 

With introducing a virtual detector as in Sec. |TT] one 
can identify eight coupling operators in the interaction 
Hamiltonian 



Ai di (g) , 
A2 = d\<E)b, 

t 



1* ,t 



A3 
A5 

Ae 
At 
Ah 



di (g) 1 , 



- ^ tlfiCkR , 
k 

k 



kRCkR • 



kL'^kL 1 



-Be - ^ tkL^^kL , 



BT = Y.t 



2* J 
feL^feL ' 



4 ® 1 , Bs^Y. ^kLCkL . 



(30) 



With introducing the continuum tunnchng rates 

k 

fe 

7fl/L(^) = '^'^J2*k*R/L4R/LH^ - <^kR/L) (31) 

we outline the derivation of the explicit Liouvillc super- 
operator in appendix IB] Note that now some tunneling 
rates may assume complex values. In the flat band limit, 
we assume them as frequency independent, which implies 
l7Lp = ^liTl2 and \jr\'^ = TriTr2. 

Generally, we find that the Liouville su- 
peroperator couples the six matrix elements 







/L(u;,/3,+f) 










FIG. 2: [Color Online] Depending on the bias voltage V , par- 
ticles may travel from left to right or vice versa through the 
effective double dot system. For negligible Coulomb interac- 
tion U one has two independent channels. 



Poo, 00, Poi.oi, Pio.iOj Poi,i07 Pio.oi- However, 

in the BMS limit (r — > oo) the populations 
Poo,oo,Poi,oi,pio,io,Pii,ii wiU - for lifted energetic 
degeneracy ei ^ £2 ~ decouple from the coherences 
Poi,iOj Pio,oi- For energetic degeneracy however, this 
decoupling does not take place in the BMS limit. The 
nondegenerate BMS case will be termed rate equation 
(4x4 Liouvillc superoperator), whereas wc will refer to 
the degenerate BMS case as quantum master equation 
(6x6 Liouville superoperator). 

Thus, one can see that the BMS Liouvillian behaves 
discontinuously with respect to the Hamiltonian param- 
eters, which is of course also reflected in observables. 
We show in appendix O how one obtains the Lamb-shift 
terms (odd Fourier transform) for Lorcntzian-shaped 
bands. 



V. BMS EQUILIBRIUM RESULTS AND 
FLUCTUATION-DISSIPATION RELATIONS 

For a system in equilibrium (either obtained by choos- 
ing identical chemical potentials in both reservoirs /i^ = 
/^/j = /i or by turning off one coupling (e.g. = Q and 
^R ~ fi or Tr ~ and /i^ = m) in the serial case) we 
obtain equilibration of both temperature and chemical 
potentials with that of the reservoir. This holds for both 
models (|25p and ((29|) . More explicitly, when the reservoir 
is in the state 



PB 



-/3(Hb-mA'b) 



(32) 



the resulting BMS Liouville superoperators has the 
steady state 



PS = 



Z's 



(33) 



where iVg ~ d\di + d\^d2 is the system particle number 
operator. The equilibration of both temperature (3 and 
chemical potential /i between system and reservoir for 
these systems is a remarkable result, since to our knowl- 
edge only equilibration of temperatures has been shown 
so far—. 

In all analytically accessible cases, we can verify the 
Johnson-Nyquist (fluctuation-fissipation) relatior*i2£ii21 



^(0) 



y=o 



2 dl 



v=o 



(34) 



which requires bi-directional n-resolved master equations 
of the form ^ and appears to be a general feature of the 
BMS approximation^^. For zero bias voltage V = we 
have shown as a sanity check that the resulting current 
vanishes. 

A fluctuation-dissipation relation for nonlinear trans- 
port based on detailed balance has recently been de- 
rived in Ref,i£^: = 1 with A = e/3{T^L, Vr} = 
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e/3{V/2, —V/2} and Q = {n^, nij} for two terminals and 
symmetric bias. We have checked this relation for a single 
noninteracting level and find it satisfied within the BMS 
approximation only up to third order in V. However, 
since we only consider counting at the right junction the 
check for the interacting systems is beyond the scope of 
the present work. 



VI. TRANSPORT RESULTS: 
CONFIGURATION 



SERIAL 



A. BMS Rate Equation 



The BMS current in Coulomb-blockade but high-bias 
limit (CBHB, formally obtained by the limits lim lim 

V — >oo U — 'oo 

such that always V -^U) equates to 



jCB ^ 



(ei - £2) Tl + {2rL + Tr) r2 



(35) 



which is also found in the numerical solution for the 
stationary current, see Fig. [31 The above result differs 
from known results in the literature^^i^^i^i by a missing 
term rx,r^/4 in the denominator. This difference is of 
higher order than the validity of the perturbation the- 
ory. However, this leads to the unphysical artifact that 
in our case for ei = £2 the current becomes independent 
on the interdot hopping rate Tc, which contrasts with 
the expectation that it should always vanish for Tc 0. 
The main reason for this failure of the BMS approxima- 
tion lies in the use of the non-localized energy eigenbasis 
during the secular approximation procedure: Tunneling 
may occur into non-adjacent eigenstates even when Tc is 
small, which finally leads to the observed current. When 
the secular approximation is not performed (technically, 
by using a finite maximum coarse-graining time) , we will 
see that the current vanishes for Tc ^ as expected, see 
Sec. IVI Bl In the actual infinite bias case (formally ob- 
tained by the limit lim such that we keep finite U) we 

obtain the stationary current 



TlTr {Tl + Tr) T^ 



(ei - £2)' TlTr + {Tl + Tr)^ T^ 



(36) 



which also displays the same artifact as the Coulomb- 
blockade current ([35])) for degenerate single-particle en- 
ergies, but is - as one would expect - symmetric un- 
der exchange of L and R. This current is also found 
in the numerical solution, see Fig. [3l and differs from^i 
by a missing term TlTr{Tl + Tr)'^/A in the denomi- 
nator. The current-voltage characteristics does now ex- 
hibit several current steps, since we explicitly allow for 
a doubly occupied system. Since we have chosen sym- 
metrized single-particle energies (£2 = ~ei)i the low- 
est two transport-channels are not distinguished. Nat- 
urally, the current-voltage characteristics is no longer 



point-symmetric for different single-particle energies. In 
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FIG. 3: [Color Online] Stationary Markovian current ver- 
sus bias voltage for the serial double dot for different single- 
particle energies. Other parameters have been chosen as 
= 1.0, U = 5.0, Tl = Tr = 0.1^ = r, /3 = 10.0. For degen- 
erate single-particle energies (black solid line) the current is 
point-symmetric around the origin. This does not persist for 
non-degenerate single particle energies and becomes strongly 
pronounced for large splittings S = ei — 62- For ei = — £2 one 
observes three non-vanishing plateau currents for both posi- 
tive and negative bias voltages (given a sufficiently low tem- 
perature). The first current step enables transport by transi- 
tions from |00) — > \E-) and |00) — > and at the second 
and third current steps the additional channels \E+) — > 
and \E-) ^ are opened. It is visible that with large level 
splitting (solid green), the second plateau equilibrates with 
the third one for positive bias, whereas it approaches the first 
plateau for negative bias. This is a localization effect (for 
explanations see the text). 



addition, when the splitting is significantly larger than 
the interdot tunneling rate Tc, we see in Fig. [3| (green 
curve) that one transport channel is strongly suppressed. 
This effect can be interpreted via the localization of the 
wave functions: For \ei — £2] » Tc (large splitting), the 
eigenvectors ([^ become localized, e.g., for £1 > £2 one 
obtains that \E+) ^ |10) and \E_) |01). When they 
are above the Fermi surface of their adjacent lead, they 
may not contribute to transport, such that only one state 
{\E-^) for positive bias and for negative bias) medi- 
ates the transport. Naturally, this also leads to asymme- 
tries in the voltage-dependent zero-frequency Fano fac- 
tor, see Fig. [4|below. Note that for extremely large split- 
tings between the single particle energies £; (constant in- 
terdot tunneling rate Tc assumed) the Coulomb-blockade 
assumption may even become invalid (not shown), as the 
energy for the doubly-charged system may fall between 
the two single-charged energies (not shown). 



From Eqn. (jl3|) wc may calculate the finite-frequency 
noise for the right lead and wc obtain for the frequency- 
dependent Fano factor Fr{lu) = Sr{uj)/\I\ in the CBHB 
regime the expression 
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WTi^Tr + TRT^^T^f + (5^ + 4r2) [5^ (ri + Tl) + (r,\ + Tl) T2] + (S^ + AT^f c.^ 



(37) 



r 



where we have used the short-hand notation 5 = ei — e2 
and Fcff = 2Tl+Tr. The zero-frequency hmit is depicted 
in Fig.m where the numerical solution coincides with the 
above result in the CBHB regime. The difference between 
the above zero-frequency Fano factor and the Fano factor 
in Ref.^ is of the order 0{TlTr/T'^). It can also be de- 
duced that the zero-frequency limit of Eqn. (|37|) becomes 
super-Poissonian (indicating bunching) 
when the right-associated tunneling rate becomes sig- 
nificantly smaller than the left-associated one Tr < 
i+2T^/i5^ , which corresponds to a highly asymmetric sit- 
uation. Assuming similar single-particle energies [5 ~ 0) 



the zero-frequency version of Eqn. (j37p reproduces previ- 
ous results in the literature^. Assuming equal tunneling 
rates F^ = F^ we see that the frequency-dependent Fano 
factor in BMS approximation has two minima only if 
the splitting between the single-particle energies is larger 
than the tunnehng rate {\5\ = |ei — 62] > V — 2\Tc\ ~ 
0.49 |Tc I) and is featureless (only one minimum) other- 
wise. This is well confirmed in the numerical solution, see 
the solid lines in Fig. [5l Even when one has two minima 
in the finite-frequency noise, Fig. [6] demonstrates that 
this only happens in the CBHB limit. For the frequency- 
dependent Fano factor in the infinite bias limit we obtain 



J 



Fr{^) = ^- 



2TlTrT^ [5^ (Fi - TlTr + T] 



(52 + 4r2) [P (F2 + F|^) + 2F2r2] w2 + (J2 



r 



(38) 



where we have again used the abbreviation (5 = ei — £2 but 
this time use F = T -\- T r. The zero-frequency limit is 
well reproduced in figure 2] in the infinite-bias regime. As 
before, the difference between the above zero-frequency 
Fano factor and the Fano factor in Refi^ is of the order 
OiT lT r/T'^)- It is straightforward to show that the zero- 
frequency version of the Fano factor psp may never be- 
come super-POoissonian. In addition, for equal tunnel- 
ing rates F^ = F^ the infinite bias frequency-dependent 
Fano factor has only one minimum - regardless of the 
splitting of the single particle energies, see the dashed 
lines in Fig. [5] 

Finally, Fig. [5] demonstrates the smooth transition of 
the frequency-dependent Fano factor from the Coulomb- 
blockade regime towards the infinite-bias limit: For small 
positive bias, the frequency-dependent Fano factor be- 
comes slightly super-Poissonian for finite uj and poisso- 
nian for w = 0, compare also Fig. ID In the CBHB 
regime (solid green line), we observe a sub-Poissonian 
Fano factor with two minima and in the infinite bias 
regime (dashed green line), the two minima merge into 
one. 



B. Coarse-Graining Results 

For finite coarse-graining times, the Coarse- Graining 
Liouville superoperator does not resemble the BMS ap- 
proximation but nevertheless preserves positivity (Lind- 
blad form), compare the discussion in Sec. IIIII We have 



seen that for degenerate single-particle energies ei = £2, 
the BMS currents (r — > 00) become independent of 
the interdot tunneling rate Tc, compare equations (j35p 
and (j36[) . That is, when Tc — ^ 0, one still has a nonvanish- 
ing current, which is completely unreasonable, as in this 
limit one has two single levels coupled to different leads. 
Note however, that the exact limit Tc ~* would in this 
case lead to a different BMS quantum master equation 
with a 6 X 6 Liouvillian that does not admit any current. 
This unphysical behavior is not found for a finite coarse- 
graining time, see figure [T] For a finite coarse-graining 
time, the stationary current vanishes quadratically as the 
interdot tunneling rate Tc goes to zero, which is con- 
sistent with earlier results for the CBHB regime in the 
Iiterature2^i2^. For larger interdot tunneling rates, the 
magnitude of the stationary current quickly saturates, 
such that the only effect of a varying Tc is the increasing 
level splitting for the single-charged states |— ) and |-|-) 
(emerging additional current steps). Therefore, for large 
Tc we obtain the previous results of the BMS approxima- 
tion (compare e.g. orange curve in Fig. [7] with the black 
curve in Fig. [3]). This behavior is quite general: When 
the coarse-graining time is larger than the inverse level 
splitting (in our case between the single-charged states), 
one can show that the coarse-graining results and BMS 
results coincide. 



Regarding 
approac h^^i^^ 



the dynamical coarse-graining 
these findings demonstrate that for 
the serial double dot considered here, the coarse-graining 
time must not be sent to infinity. To obtain non- 
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e^=M).S e,=-0.5 : 
+1.0 e'=-1.0 
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bias voltage [a.u.] 

FIG. 4: [Color Online] Zero frequency Fano factor versus 
bias voltage for different single-particle energies. Parameters 
and color coding coincide with Fig. [S] For degenerate single- 
particle energies (black) we have a symmetric Fano factor. For 
large splittings 5 = £1—62 (green) the equilibration behavior of 
the Fano plateaus coincides with that of the current plateaus. 
The peak at the origin results from (thermal) Nyquist noise, 
cf.^. 
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0) 



FIG. 6; [Color Online] Frequency dependent Fano factor 
Fe.{uj) for large splitting versus bias voltage and lo. Other 
parameters have been chosen as in Fig. O but the single- 
particle energies have been fixed at e\ — +1.0 and £2 ~ —1.0. 
In the CBHB regime (e.g. V — ±5), there are two valleys 
corresponding to two minima of the Fano factor, whereas this 
feature is lost in the infinite bias regime (e.g. V = ±20), 
compare also Fig. O 
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FIG. 5; [Color Online] Frequency dependent Fano factor 
Fft,{u) for different splittings in Coulomb-blockade high-bias 
(CBHB, solid) and infinite bias (dashed) limit versus u. 
Other parameters have been chosen as Tc = 1.0, U = 5.0, 
Vl = Tr = 0.1 = r, and (3 = 10.0. The frequency-dependent 
Fano factor exhibits two minima that are for large splittings 
5 S> Tc situated at uj* = ztF (solid lines). In the infinite bias 
limit, this feature is generally lost (dashed lines). 



Markovian effects, it must for small times scale linearly 
with the physical time t = t, but for large times it must 
saturate r = Tmax in order to avoid artifacts such as 
currents through disconnected structures. 



VII. TRANSPORT RESULTS: PARALLEL 
CONFIGURATION 

A. BMS Rate Equation 

When the single-particle energies are non-degenerate 
£1 7^ £2, the equations for the populations decouple from 
the equations for the coherences, see appendix [B] for de- 
tails. For these rate equations we obtain in the CBHB 
limit the current 



CB 



(rii + Tl2) ri^irij;2 



R2 



Tt.oT 



L2i Rl 



flli R2 



(39) 



For equal left- and right-associated tunneling rates 
(Fii = Tl2 = Tl, Tri = Tr2 = Tfl) we recover the 
known resulO^a 2rLri^/(2rL + Ffl). Clearly, the above 
current vanishes if just one of the right-associated tun- 
neling rates goes to zero. This is a Coulomb-blockade 
effect: Due to the high-bias limit, an electron gets stuck 
in the orbital with the vanishing right-associated tunnel- 
ing rate, compare also Fig. [21 This electron blocks the 
transport through the other channel by the Coulomb in- 
teraction. The current ((39)) is also well confirmed by the 
numerical solution, see Fig.[8l In contrast, in the actual 
infinite bias case we obtain the current 



rLirL2 (r_Ri -i- ri^2) + {Tli + ^l2) r_Riri^2 

{TLi+Tm){TL2 + TR2) 



,(40) 



which is again completely symmetric under exchange R 
and L. When one left- or right-associated tunneling 
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FIG. 7: [Color Online] Plot of the stationary current versus 
bias voltage for different interdot tunneling rates Tc. The 
coarse-graining time has been chosen as r = 100, the inverse 
of the coupling strengths Fl and F/j. For small Tc, we observe 
a quadratic growth of the current magnitude (dashed curves) , 
which quickly saturates. For larger Tc (solid lines), the current 
magnitude becomes independent of Tc as also seen in the BMS 
approximation for ei = £2, but the level splitting between the 
single-charged states |— ) and |-|-) is slowly resolved (emerging 
additional current step). Other parameters were chosen as 
ei = £2 = 0, i7 = 5.0, Fi = Ffl = 0.1^ and P = 10.0. 



rate vanishes, we obtain the current of a single-resonant 
level^^, which is independent of the remaining tunneling 
rate associated with the blocked channel: As the infinite 
bias limit explicitly allows double occupancy, we do not 
observe Coulomb blockade in this case. Also the infinite- 
bias current is well reflected in the numerical solution in 
Fig. [HI Wc sec that for the rate equations, the allowance 
of doubly occupied states simply opens more transport 
channels, which are visible by at most four steps in the 
current-voltage characteristics. Interestingly, the height 
of these steps is independent on the level splitting, but 
is just determined by the tunneling rates. However, the 
(sequential tunneling) excitation spectrum of the system 
Hamiltonian can be probed by the position of the cur- 
rent steps (sufficiently low temperatures provided) . Note 
that the current for Si = £2 provides the analytic contin- 
uation for nearly degenerate energies. In the degenerate 
case the quantum master equation with coherences has 
to be applied as it will be shown in Sec. I VII 5] However, 
when additional dephasing processes lead to a decay of 
coherences, the rate equation approach can still be mo- 
tivated. 

For the frequency-dependent Fano factor we obtain in 
the Coulomb-blockade high-bias limit 
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FIG. 8: [Color Online] Stationary Markovian rate equation 
current versus bias voltage for the parallel double dot for dif- 
ferent single-particle energies. Other parameters have been 
chosen as U = 5.0, Fii = = Fm = 1^r2 = 0.1^ = F, 
/3 — 10.0. For sufficiently low temperature, the current steps 
may be found at the (doubled) excitation spectrum of the 
system Hamiltonian (sequential tunneling limit). As all tun- 
neling rates are assumed equal, the curves are always point- 
symmetric around the origin. For sufficiently large splittings 
(green curve) and low temperatures one may resolve between 
different transport channels. The height of the plateaus is 
independent of the system energies. For symmetric single- 
particle energies one has three plateaus with a non-vanishing 
current. With increasing splitting (green), the second plateau 
broadens thereby annihilating the first one. This happens 
when the state |11) becomes energetically more favorable than 
1 10), such that the coulomb-blockade assumption is not ap- 
plicable. 



where we have chosen the abbreviations 



Mo = 
Vi = 



^L2 (Zr^i -|- Tl2)^bj — Ll^ L2^ Rl^ B2 

+ (rii + 2riiri2 + r|ji)rL, 

rLi' + r^i' + 2rii (ria + r^i) 

-(rL2 + rH2)'l . (42) 



(41) 



For zero frequency and similar left and right tunneling 
rates (F^i = Tl2 and Tjn = Tfi2) this coincides with 
previous resultsS^, see also table I(column (ii)) in RefiiS^. 
As expected from the current in the CBHB regime, it can 
be checked easily that the zero-frequency version of the 
above Fano factor is always super-Poissonian when one 
of the right-associated tunneling rates vanishes (bunch- 
ing due to dynamical channel blockade) and is always 
sub-Poissonian when one of the left-associated tunneling 
rates vanishes. Super-Poissonian Fano factors for asym- 
metric systems are well-known in the literaturoSS,. The 
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FIG. 9: [Color Online] Zero frequency Fano factor (right) 
versus bias voltage for the parallel double dot for different 
single-particle energies. The parameters correspond to those 
of Fig. [S] The peak of the zero-frequency Fano factor at small 
bias voltages results from thermal (Nyquist) noise. 



corresponding infinite-bias Fano factor equates to 



F^{lu) = 1 - 2 X 



(43) 



F Ll^ Rli^ L2+T R2) 



(ri,2+r^ 



For zero frequency and equal left and right tunneling 
rates (Fli ^ Tl2 and Tm = Tii2) this also coincides 
with previous results^^, see also table I (column (iv)) 
in RefiiSS. From the positivity of the tunneling rates 
it is obvious that the above Fano factor is always sub- 
Poissonian, which well matches our expectations for the 
infinite bias limit. These results are well confirmed by the 
numerical solution, see Fig.[n]for the zero-frequency Fano 
factor. Note that we consider only the right-associated 
noise here, which only at zero frequency coincides^^ with 
the full frequency-resolved noise. 



It follows from equations (^1]) and ([33]) that for equal 
tunneling rates the frequency-dependent noise in the 
high-bias coulomb blockade and in the infinite bias limit 
will essentially be featureless with a single minimum. Nu- 
merically, we see that for equal tunneling rates this ex- 
tends to the whole bias range (not shown). However, 
when we assume varying tunneling rates, we may find 
features in the frequency-dependent noise for the CBHB 
regime - this can be found analytically from Eqn. (|¥T|) . 
These findings suggest that a feature in the frequency- 
dependent noise is not a unique indicator for coherent 
dynamics in the system. This would require the evalua- 
tion of the full noise, which is the weighted sum of noise 
at left and right junction and cross-correlators, see e.g. 
in Ref42.. This is beyond the scope of this paper. 



B. BMS Quantum Master Equation 

When the master equation is considered, we assume 
energetic degeneracy ei = €2 = e throughout. For sim- 
plicity, we assume 7l = V^liTl2 and -jr = VTmHu, 
which corresponds to the fiat band limit. In the infinite- 
bias limit, we obtain the general stationary master equa- 
tion current 



r 



^Ll^Rl + ^L2^R2 + 2\/T LiT Rl^ R2 



LI 



L2 



Rl 



+ T 



R2 



(44) 



which does not necessarily coincide with the correspond- 
ing infinite bias rate equation result (j40[) . For gen- 
eral tunneling rates we may also evaluate the frequency- 
dependent Fano factor analytically, which is however too 
lengthy to be reproduced here. However, when we as- 
sume some symmetries, the expressions simplify signif- 
icantly. For example, when of the four tunneling rates, 
two pairs are mutually equal, one has three different com- 
binations, which we will label as (compare also Fig. ^ 

I. left-right symmetric coupling Fli = Tri = Fi with 
^L2 = ^R2 = r2, 

II. top-down symmetric coupling F^i = Tl2 = 
with Tm = Tr2 = r_R, and 

HI. antisymmetric coupling F^i = Tfj2 = with 
Tl2 ^ ^Ri = Tb (compare also ''^i^° ) 

in the following. Note that naturally, when all tunneling 
rates are the same, the above three cases coincide. In 
addition, we note that in all three cases, the infinite bias 
master equation current (j44p and the infinite bias rate 
equation current (|40)) coincide. 



Symmetric configurations: I and II 

Under the symmetry assumptions I and II, the 
frequency-dependent infinite bias Fano factor reduces to 



2{Ti+r2f+u;^ 



4(ri 
1 - 



4 (Fl + Ffl)' + 



(45) 



Also when the corresponding symmetries are used, these 
frequency-dependent Fano factors differ at finite fre- 
quency from the corresponding rate-equation result (|43p . 
the Fano factors only coincide at zero frequency. 

It should be noted that for exact symmetry as- 
sumptions I and II, the null-space of the Liouville 
superoperator L{Q) becomes two-dimensional for the 
complete bias range. (Consequently, the steady state 
and the observables would depend on the choice of 
the initial state in that case.) However, when we 
also consider couplings between the subspaces, we 
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obtain an unique steady state, which is for symme- 
try assumptions I and II analytically continued into 
Eqns. (j45|) . When symmetry assumptions I and II coin- 
cide (all tunneling rates equal), this becomes obvious in 
the basis poo,oo, 5 [poi,oi + Pio,io + Poi,io + Pio.oi], 

5 [POl.Ol + PlO.lO ^ POl.lO ^ PlCOl]; 

^ [poi,oi - Pio,io]: ^ [poi,io - Pio,oi], where the 
Liouvillc superoperator L{x) decouples into three 2x2 
blocks, two of which have a null-space as x ^ 0- 
Therefore, the above infinite bias Fano factors should be 
interpreted as being valid for approximate symmetries I 
and II. For exactly fulfilled symmetries one will obtain 
two partial currents and Fano factors. Since the finite- 
bias behavior is especially interesting, we will discuss 
these subtleties in the following at 1/ = [/ : In addition, 
we assume vanishing single-particle energies ei = 62 = 
for simplicity. When one is faced with a two-dimensional 
null-space of the Liouville superoperator, one has to 
choose the null-space vectors orthogonal and of course 
normalize with respect to the trace. Then, we obtain 
the currents in the two subspaces (A) and (B) 
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for symmetry assumption I and similarly for symmetry 
assumption II 
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Naturally, for completely symmetric tunneling rates, the 
currents in the respective subspaces (A) and (B) of ([46|) 
and (|T7)l coincide. Regarding the noise, we obtain the 
frequency-dependent Fano factors for symmetry assump- 
tion I 
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Vi ^ {y - 1) (1 + y^) with y = ef^^^ in the last line. 
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where we have used in the first equation Co 



[4(F 



'] and in the last line A/q 



4ri [2rLTRro + (Fi + F|) Vi] , W = 2rRy^ + TlVi 
with the polynomials Vq, = l-2y + Ay'^ - 2y^ + y*,Vi = 
1 — y + 2y^ — y"^ + y'^, and also as before y = e^^^"^. Also 
for the Fano factors, the expressions in equations ([48|) 
and coincide in their respective subspaces (A) and 
(B) for completely homogeneous tunneling rates. 



Antisymmetric configuration III 

Under the assumption of an antisymmetric configura- 
tion we obtain an unique steady state - with the excep- 
tion of completely homogeneous couplings. 

The infinite-bias Fano factor reads 
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-Tb)', Vi 
(F2+6F,Ffc + F2), 
coincides with the 
43]) only at zero 



corresponding rate equation result 
frequency. 

The situation becomes sophisticated when we consider 
the finite bias case V = U together with symmetry as- 
sumption III, where (Fq 7^ F^ provided) one only has 
one stationary state - which is also the generic case for 
arbitrary couplings. In this case, the Lamb-shift terms 
directly affect the current: Formally, this is visible in 
the appearance of Digamma functions '^(x) (resulting 
from the Cauchy principal value integrations, see ap- 
pendix [C|) in the stationary current. With the replace- 
ments y = f = Ta + Tb, f] = Ta ~Tb and 
^{y) = [* (i + zii^) - * (5 - 3*^)] we obtain 
with symmetry assumption III ai V ~ U and ei = £2 = 
for the current after a tedious and lengthy calculation 
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bias voltage [a.u.] 

FIG. 10: [Color Online] Markovian stationary current pre- 
dicted by the BMS quantum master equation for exact de- 
generacy ei = 62 = e. The parameters that are the same for 
aU curves are U = 5.0, Tli = 0.1^ Tm = 0.15^ Vr2 = 0.1^ 
j3 = 10. The other parameters have been varied as follows: 
Case A (black): e = 0.0, Tl2 = 0.15^ Case B (red): e = 1.0, 
Vl2 = 0.15^ and Case C (blue): e = 0.0, V l2 = 0.05^ The 
thin dashed lines show the corresponding rate equation re- 
sults. The nearly complete suppression of the current for 
large PU at 1^ « (7 is for exact degeneracy quite robust with 
respect to the remaining parameters. 



where we have abbreviated the polynomials Vi = 1 + y 
and V2 = 1 ^ ?/ + y^- This is consistent with previous 
numerical results^. For low temperatures (large /3J7), 
the above expression behaves asymptotically as 



rv=u PUJ^^ SraFt 



(52) 



i.e., the current is exponentially suppressed for low tem- 
peratures! The numerical solution confirms this result 
and shows that it goes along with a negative differential 
conductanc e^^ '^'^ , which is - exact energetic degeneracy 
ei = £2 provided - quite robust with respect to the re- 
maining parameters, see Fig. [TOl For example, assuming 
non-vanishing single-particle energies ei = £2 7^ 0, we 
see that the current suppression valley is simply shifted 
away from the origin. Even when the symmetry assump- 
tions regarding the coupling strengths are not obeyed, 
the current suppression is qualitatively robust. When 
we evaluate the zero-frequency Fano factor, we observe 
huge Fano factors in the current suppression region, see 
Fig. [TT] Note however, that super-Poissonian (but sig- 
nificantly smaller) Fano factors may also be observed in 
the rate equation case (due to dynamical channel block- 




bias voltage [a.u.] 

FIG. 11: [Color Online] Zero frequency Fano factor pre- 
dicted by the quantum master equation for exact degeneracy 
ei = £2 = e. Parameters and color coding are the same eis 
in Fig. 1101 The current suppression in the quantum master 
equation goes along with highly super-Poissonian Fano fac- 
tors. 



ade as discussed before), compare also the dashed lines 
in Fig. [ini 

At the point of completely homogeneous tunneling 
rates, we can obtain currents and Fano factors for the 
whole bias range 
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where we have abbreviated x = e^^l"^ and y = 1"^ in 
the last line, see Fig. [T^l For V — U , the frequency- 
dependent Fano factors can be derived from (j48p when 
Fi = Fz = F or from (gSl) when Vl=^r = F. 

The alert reader will have noticed that in case of com- 
plete coupling symmetry F^i = Vi^2 = ^ri = ^R2 = T 
the currents // and /// and the Fano factors Frj{uj) and 
Frji{uj) coincide in their respective subspaces, compare 
equations and (|T7)) . However, the antisym- 

metric current lY/j^ does neither converge to the current 



15 



CD 

c 

,0 -0.005 





1 




1 




: r=o.r 








i„:r=o.r 








r^=o.i5 




1 it 


- 


r,-o.ir 




i / ' 




















/" 
/ 

1 / 




- 




If! 

if: 


■■J 






if: 


/ 




















1 




1 



bias voltage [a.u.] 

FIG. 12: [Color Online] Plot of the stationary currents I a and 
Ib for the decoupling subspaces for F = 0.1^ together with the 
(unique) stationary current obtained for Fa = F = 0.1^ and 
different values of Fj, versus bias voltage. Other parameters 
have been chosen as (3 — 10.0, U — 5.0, and ei = £2 = 0. For 
nearly equal tunneling matrix elements, the actual current 
(solid line) is a convex superposition of the partial currents 
(dashed lines). 
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FIG. 13: [Color Online] Logarithmic Plot of the Fano fac- 
tor Fr{0) versus bias vohage -20.0 <V< +20.0 for differ- 
ent values of F;,. Parameters and color coding correspond to 
Fig. 1121 When the partial currents do not coincide, the Fano 
factor diverges for Fa = Ft, i.e., it cannot be written as a 
convex combination of the partial Fano factors (dashed red 
and blue lines). 



of subspace (A) nor to the current of subspace (B) when 
complete symmetry is assumed, compare Eqn. (|5ip . This 
demonstrates that for near degenerate tunncUng rates, 
the actual current will rather be a superposition of the 
partial currents. Indeed, one can find convex linear com- 
binations of the two stationary states that reproduce the 
current ////jp p. compare also Fig. 1121 An interest- 
ing consequence is that in this case (decoupling subspaces 
with different partial currents), it follows directly from 
the definition of the cumulant generating function ([7]) 
that the (zero-frequency) Fano factor nmst diverge (tele- 
graph noise), which is also seen in the numerical investi- 
gations, see Fig. [T31 Note that this divergence will also 
hold for the higher cumulants^. 

It should be stressed that the huge Fano factors ob- 
served in Fig. [T3] around the Coulomb-blockade region 
are obtained for highly symmetric systems. For reason- 
able parameter values, they are by orders of magnitude 
larger than the super-Poissonian Fano factors observed 
for asymmetric systems^^. In the frequency-dependent 
Fano factor, this divergence goes as a very narrow and 
tall peak within a large valley at low frequencies, see 
Fig. [T31 This is similar to the spectral Dicke effec t^^i^^" . 
An experiment will average over a finite frequency inter- 
val determined by the inverse measurement time. Since 
the area below the peak remains approximately constant 
(compare the scaling of height and width with respect 
to the asymmetry rj in Fig. I14p . the super-Poissonian 
behavior can in principle be resolved. However, when 
the frequency interval of the measurement is too large, 
the valley may dominate the super-Poissonian peak and 
one may obtain a sub-poissonian Fano factor. Therefore, 
the inverse measurement time should scale with the peak 
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FIG. 14: [Color Online] Frequency-dependent Fano factor for 
symmetry assumption III at the maximum current suppres- 
sion V — U for different coupling asymmetries rj. Parameters 
have been chosen as ei = £2 = 0, /? = 10, C/ = 5, Fa = F — rj, 
Fi, = F -I- 77, where F = 0.1^. The divergence of the actual 
zero-frequency Fano factor for 77 — > manifests itself as an 
extremely slim peak within a large valley around the origin. 
This is similar to the spectral Dicke effectiii^. 



width. 



C. Telegraph statistics - Distribution function P„{t) 

It is quite instructive to study the impact of these 
results on the time-dependent probability distribution 
Pn(t). Given the cumulant-gencrating function, it is un- 
fortunately non-trivial to obtain P„(i) via the inverse 



16 



Fourier transform of equation since that would in- 
volve an integral over a highly oscillatory function. How- 
ever, when the cumulant-generating function is highly 
peaked in the interval x S (— tt, -I-tt), we may calcu- 
late it approximately via the saddle-point approxima- 
tion (also termed stationary phase approximation, sec 
e.g. Refs^^^iii) 



+7r 



2^ 



,S"(x*,t)(x-X*)V2 



dx, (54) 
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where we see that the remaining integral in the above 
equation just corresponds to a normalization of the dis- 
tribution, since it docs not depend on n. The position 
X* of the integrand peak is determined by the equation 

d^S{x,t\^^,=in, (55) 

which only admits purely imaginary solutions for x* ■ 
Since Eqn. (j55p is rather demanding to solve numerically, 
we have computed n from Eqn. (j55p and the correspond- 
ing distribution Pn{t) ex. cxp{S'(x*,t) — inx*} paramet- 
rically as a function of all imaginary x* ■ Normalization 
was then performed afterwards, see Fig. 1151 Instead of 
obtaining a bimodal distribution as might be naively ex- 
pected from the existence of two different currents in the 
subspaces A and B, we see that the actual probability dis- 
tribution is unimodal^. The mean of the distributions 
- divided by the snapshot time t - yields the current at 
the corresponding bias voltage, which is well consistent 
with the current-voltage characteristics in Fig. [121 The 
huge zero-frequency Fano factors for bias voltages below 
the Coulomb interaction strength - compare Fig.[T3]" are 
reflected in extremely large widths for the corresponding 
bias range. Note however that Fig. [T51 shows a finite time 
snapshot which does not capture the large time limit of 
the zero-frequency Fano factor. 

D. Coarse-Graining Results 

Although the BMS currents do not display completely 
unphysical behavior, it is disturbing that the BMS Li- 
ouvillian behaves discontinuously as a function of its pa- 
rameters, which also transfers to observables such as the 
current. 

One may argue that in a realistic setting, there will 
always exist lifted degeneracies, such that for the sta- 
tionary state the result of the rate equation is relevant. 
In this case, the DCG approach will finally approach the 
BMS rate equation current, but the BMS master equa- 
tion current will appear as a metastable state - regard- 
less of the initial state. These expectations are also found 
in the numerical solution, see Fig. 1161 For small times 



FIG. 15: [Color Online] Normalized probability distributions 
P„{t) for different bias voltages versus the number of tun- 
neled particles n obtained via the saddle-point approxima- 
tion (|54p . Other parameters have been chosen as [/ = 5, 
/3 = 10, t = 20000, and r„ = r - 77, Ft = r + 77, where 
r = 0.1^ and r; — 10~'' . The inset relates the probability dis- 
tribution for V — 9.7 with the corresponding partial currents 
(lines with symbols) that one would obtain in subspaces A 
and B for ri — (cf.— ). For bias voltages near the current 
suppression, the distribution becomes extremely flat - consis- 
tent with the huge super-Poissonian Fano factors. In contrast, 
for large bias voltages we approach the infinite bias distribu- 
tion, which is characterized by a finite width. For V = U, 
the exponential suppression of the current is manifest as an 
exponentially decaying distribution (straight line). 



t = T, the transient behavior resulting from the initial 
state is still found in the stationary current (solid lines 
in Fig.[Tn]), but the relaxation drags the current towards 
the master equation result. For times larger than the in- 
verse level splitting, we observe a relaxation back towards 
the rate equation result (which had also been chosen as 
the initial state in Fig. \W\\ . However, as the timescales 
of relaxation (F^^ = 100) and the inverse level splitting 
[S^^ = 1000) do not completely separate, we do not ob- 
serve a complete decay into the master equation result. 
The inset in Fig. 1161 shows that this behavior is indepen- 
dent of the particular initial state and that for separating 
time scales a complete decay into the metastable state is 
observed. 

Inspired by the success of finite maximum coarse- 
graining times ((2T|) for the serial model in Sec. IVIBl we 
have also calculated the stationary current as a function 
of the coarse graining time. We find that the finite coarse- 
graining (FCG) approach also predicts the qualitative ef- 
fects of the BMS master equation in the stationary state, 
see Fig. 1171 By varying the coarse-graining time we actu- 
ally observe a smooth crossover from the quantum mas- 
ter equation results (small coarse-graining times) towards 
the rate equation results (large coarse-graining times). 
This directly demonstrates that the singular coupling 
limit is not the only master equation method yielding 
Lindblad type master equations that smoothly interpo- 
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FIG. 16: [Color Online] Plot of the time-dependent DCG cur- 
rent versus bias voltage for different times. The other param- 
eters were chosen as ei = -1-0.0005, £2 = —0.0005, U = 5.0, 
(3 = 10.0, Tli = Tfla = 0.1^ and r^a = Tm = 0.15^ In 
the plot, the initial condition has been chosen as the sta- 
tionary density matrix of the BMS rate equation, whereas 
the inset shows the actual time-dependence of the current at 
V — U (compare the vertical green dashed line) for different 
initial conditions: the steady state density matrix of the BMS 
rate equation (SS-RE), the stationary density matrix of the 
BMS master equation (SS-QME), the empty system, and the 
doubly occupied system. For times smaller than the inverse 
level splitting AE~^ = 1000 (thin solid lines), we observe a 
transient crossover from the BMS rate equation current (bold 
dashed red line) towards the BMS master equation current 
(bold dashed blue line). When the time approaches and ex- 
ceeds the inverse level splitting, this trend is reversed (thin 
dashed lines), until one recovers the initial dynamics. The 
inset shows that this behavior is qualitatively independent on 
the initial state (IS) and that for significantly smaller split- 
tings = 10^ (solid lines) in comparison to 5^^ = 10"^ 
(dotted lines) the two transient processes do well separate. 



lates between rate equation and quantum master equa- 
tion results^i. Note that alternatively, we could have 
taken the coarse-graining time as constant and modified 
the level splitting £1 — 62, compare also Fig. [71 

Provided that the measurement time exceeds the (non- 
vanishing) inverse level splitting, experiments should be 
able to back up cither the FCG results (e.g., when a neg- 
ative differential conductance is measured) or the DCG 
results (steplike current voltage characteristics) just by 
measuring the current in a fixed time interval. 



VIII. SUMMARY 

We have studied single-electron transport through se- 
rial and parallel configurations of two interacting levels. 
We considered the weak-coupling limit with arbitrary 
Coulomb interaction strengths and bias voltages. Our 
calculations in the energy cigenbasis included couplings 
between populations and coherences, which are mediated 



FIG. 17: [Color Online] Plot of the stationary current versus 
bias voltage for different coarse-graining times r. For small 
coarse-graining times, we recover the qualitative behavior of 
the BMS master equation, whereas for large coarse-graining 
times, the BMS rate equation result is approached. The other 
parameters have been chosen as ei = -1-0.0005, £2 = —0.0005, 
U = 5.0, Tli = Tr2 = 0.1^ Pia = Pm = 0.15^ and (3 = 
10.0. 

by the Lamb-shift. 

We have derived the n-resolved quantum master equa- 
tion by virtue of an auxiliary detector at one junction. 
Obtaining the Liouvillc superoperator by coarse-graining 
generally leads to Lindblad-type master equations. The 
BMS approximation naturally arises in the limit of infi- 
nite coarse-graining times. As a remarkable property, we 
find that it leads to equilibration of both temperatures 
and chemical potentials between system and reservoir for 
zero bias. Unfortunately, the BMS Liouvillian depends 
discontinuously on the parameters of the system Hamil- 
tonian: For a degenerate spectrum, coherences couple to 
the populations in the system energy eigenbasis, whereas 
for a non-degenerate system spectrum, coherences and 
populations evolve independently. We demonstrate that 
for finite coarse-graining times this is not the case. 

In the serial configuration, we observe a steplikc in- 
crease of the current versus bias voltage caused by addi- 
tional transport channels beyond the Coulomb-blockade 
regime. However, the neglect of coherences is only valid 
for large values of the interdot tunnel coupling: For van- 
ishing intcrdot tunneling, the BMS current remains fi- 
nite, which is clearly unphysical. With coherences (finite 
coarse-graining times), the stationary current vanishes 
quadratically for small interdot tunneling in agreement 
with exact results. In the Coulomb-blockade regime, the 
noise spectrum at one junction may display additional 
minima, and the zero-frequency Fano factor indicates 
bunching for highly asymmetric lead couplings. 

In the parallel case, one may also obtain super- 
Poissonian Fano factors for the rate equation case (ne- 
glect of coherences) for asymmetric lead couplings due 
to dynamical channel blockade. In contrast, the quan- 
tum master equation predicts giant Fano factors for sym- 
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metric configurations due to interference effects. This 
goes along with exponential current suppression for 
sufficiently low temperatures in the Coulomb-blockade 
regime. For perfectly symmetric lead couplings, the 
quantum master equation dynamics decouples into two 
independent subspaces of physical relevance, which bears 
strong similarities to classical telegraph signals. In the 
frequency-dependent noise spectrum, this zero-frequency 
divergence appears as a i5-like peak, similar to the spec- 
tral Dicke effect. The time-dependent DCG current dis- 
plays intriguing non-Markovian features such as the tem- 
porary dwell in a metastable state before eventually re- 
laxing into the BMS steady state. However, the final 
decay will only take place when the maximum coarse- 
graining time is larger than the inverse splitting between 
the singly-charged cigcnstatcs. 

Future research should clarify the question which max- 
imum optimal coarse-graining time should be chosen in 
order to obtain the best stationary observables. I.e., so 
far we only have discussed the effect of finite coarse- 
graining times qualitatively. A direct comparison of the 
FCG approach within exactly solvable models might give 
a quantitative estimate on the maximum coarse graining 
time Tniax and the evolution T{t) [sec Eq. (Hi}]. 
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APPENDIX A: SERIAL COARSE-GRAINING 
LIOUVILLIAN 

By transforming the bath coupling operators in equa- 
tions (|27[) to the interaction picture we obtain with using 
the tunneling rates ((28|) the bath correlation functions. 
The Fourier transforms (jl9p of the nonvanishing bath 
correlation functions arc simply given by 

712(0;) = TL(uj)fL{(^) , 

721M = TLi-Lu)[l- hi-i^)] , 

734(0)) = TR{uj)fR{uj) , 

743(0^) = Tr{-u;)[1- fRi^Lu)] . (Al) 

When we make use of the relations 
0(^1 — ^2) = I [1 + sgn(ii — ^2)] and also 

0(^2 — ii) = 5 [1 — sgn(fi — ^2)], we can insert the 
even and odd Fourier transforms ^9)) into the coarse- 
graining Liouvillian (|17p . The Liouvillian separates into 
a dissipative (oc jij{i-u) and a unitary (Lamb-shift) part 
(cx CTy(o;)), 
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In appendix [C] we demonstrate how for Lorentzian- e.g. as 
shaped bands r(o;) the odd Fourier transforms 
can be extracted analytically from the even Fourier trans- 
forms 7y(o') given in Eqn. jATJ. The time-dependence of 
the system operators arises from the interaction picture, 
it can always be written as a sum over oscillatory terms 

(eigenoperator decompositioiiS2*8i). Therefore, the time and similarly for the other integral, where the band filter 
integrations in Eqn. (|J2|) can be performed analytically, function sinc(a;) = has been introduced. Due to 

the two time integrations, products of two band filter 



19 



functions arise, and in the large r limit we may use for 
discrete 0Ja,0Ji, the identit}»2i 



T r T ' 

A(wa,t^6) = lim —sine {uj - uja)- 



smc 



5{UJ - LUa) : 



(A4) 



where Si^^^^^^ denotes the Kronecker symbol and 5{Lu — LUa) 
the Dirac S distribution. Thus, in the limit t oo also 
the frequency integral in (jA2| collapses and we obtain 
the BMS Liouvillian. 



As these Fourier transforms directly determine the 
dissipative part of the Liouvillian, this also demon- 
strates that one cannot always identify the latter 
with the real part of the Liouvillian (the tunnel- 
ing rates 7fl(w) and 7l(w) may be chosen complex- 
valued). Making use of 6(ii — = + sgn(ti — t2)] 
and Q{ti — t2) = 5 [1 ~ sgn(ti — t2)] we can insert both 
even and odd Fourier transforms (jl9p into the coarse- 
graining Liouvillian p7p . This leads to 



APPENDIX B: PARALLEL COARSE-GRAINING 
LIOUVILLIAN 



To obtain the interaction picture, it is advantageous^ 
to expand the system coupling operators in ([30]) into 
eigenoperators of the system Hamiltonian, as for ex- 



ample di (t) ~ 

e-*('^i+'^)*di4d2 
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g-ieit^^^^^^ and similarly for the 
other operators. With the tunneling rates ((3T|) wc ob- 
tain for the Fourier-transforms ([TO]) of the non- vanishing 
bath correlation functions the result 
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We demonstrate in appendix [C] how for Lorentzian- 
shaped bands r(w) the odd Fourier transforms (7y (w) 
can be extracted analytically from the even Fourier trans- 
forms 7y(tf) given in Eqn. (|Bip . Similarly to the discus- 
sion in appendix [X] we may obtain the BMS limit ana- 
lytically by letting t 00. 



APPENDIX C: EVALUATION OF LAMB-SHIFT 
TERMS 



Given the Fourier transform of the bath correlation 
functions 7^ (w) , we can obtain the odd Fourier transform 
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for the conventions chosen in (1191) via 



+00 



+00 



+00 

TT I UJ — U! 



duj , (CI) 




where P denotes the Cauchy Principal Value. Thus, 
the Lamb shift is equivalent to the exchange field in 
Refs^i23. With the relation tanh(7r?/) = 2/7r3*(l/2 + 
iy) (where "^{z) denotes the Digamma function) we can 
write 
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FIG. 18: [Color Online] Poles of the function x[^')/W - ^) 
with x(i^') = |[*+|*(| - i '^'"^"^' )] f (o)') in the complex 
plane. The Digamma function contributes poles on the lower 
imaginary axis, whereas the Lorentz function (|C5|I has two 
complex conjugate poles of first order, one of which lies within 
the integration contour. The integral along the large half 
circle vanishes due to the Lorentzian cutoff, the integral along 
the real axis corresponds to the Cauchy Principal value of 
the left hand side of Eqn. HC4p . and the integral along the 
small half circle yields the first term on the right hand side 
of Eqn. HC4|I . The situation can be treated analogously for 

^X{^)^[l-f{~L0)]f{-L0). 



which is valid for all u) on the real axis. 

Since the integral in Eqn. (jCip goes along the real axis 
only, we can exploit the above identities in the Kramers- 
Kronig relation 
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-27ri Res — 



X(t^') 



which is valid whenever the (holomorphic) function x('^) 
decays sufficiently fast in the upper complex half plane 
and where ujk denote the poles of x{^) iii the upper com- 
plex half plane. Now, the imaginary part of the above 
Kramers-Kronig relation reads 
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UJ — UJ 



(C4) 



+27r5R y Res^ 



U)'=LjJk 



and can be directly related to Eqn. (jCip when we choose 



the function x(ijj') = 
see also Eqn. ()C2 



27r 



)] f (±C.'), 



The additional factor V{±lo') 
arises from the (in reality) frequency-dependent tun- 
neling rates, compare the Fourier transforms in equa- 
tions (|Aip and (|Bip . respectively. For simplicity, we take 
Ta{Lo) = Ta^{uj) and assume a Lorentzian dependence 



The Digamma function '^{z) has poles for non-positive 
integers. Therefore, when "^ui > (the upper complex 



half plane) the function 



2tt J 



de- 



fined in Eqn. (jC2[) would behave analytically, such that 
one would not have to evaluate residues. The Digamma 
function alone would not decay asymptotically, but con- 
vergence is ensured by the Lorentzian cutoff \Cb\ . This 
Lorentzian cutoff leads to a single pole aX uji = e + i5 in 
the upper complex half plane. In principle, more compli- 
cated spectral densities can be fitted by a sum of many 
LorcntziansSI. The situation is depicted in Fig. [TS] 

With these considerations we can express the Lamb- 
shift terms in Digamma functions, where the Fourier 
transforms in equations (jAip and (jBip lead to terms like 



+00 



+00 I 
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and also terms like 

+00 



+00 I 
2 



IT \ Z ZTT 



-did 



6 + i{e + w) 



For nearly flat tunneling rates we observe a logarithmic 
divergence as 5 — > 00. These logarithmic counterterms 
can not always be neglected: For the parallel configu- 
ration, we do indeed observe their cancellation in the 
Liouvillian, whereas for the serial configuration they do 
not cancel. In the latter case, the neglect of these terms 
would even lead to negative density matrices. We suspect 
that the reason for this lies in the discrepancy between 
the localized and the energy-eigenbasis. 
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